Formation of Electronic Nematic Phase in Interacting Systems 
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We study the formation of an electronic nematic phase characterized by a broken point-group 
symmetry in interacting fermion systems within the weak coupling theory. As a function of inter- 
action strength and chemical potential, the phase transition between the isotropic Fermi liquid and 
nematic phase is first order at zero temperature and becomes second order at a finite temperature. 
The transition is present for all typical, including quasi-2D, electronic dispersions on the square 
lattice and takes place for arbitrarily small interaction when at van Hove filling, thus suppressing 
the Lifshitz transition. In connection with the formation of the nematic phase, we discuss the origin 
of the first order transition and competition with other broken symmetry states. 
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I. INTRODUCTION 

Recently, there have been reports of experimental ev- 
idence of inhomogeneous and/or anisotropic quantum 
ground states in strongly correlated systemsiSii^ The 
inhomogeneous phase — dubbed the stripe phase — was 
observed in high temperature cuprates via elastic neu- 
tron scattering experiments & Theoretical studies of the 
inhomogeneous and anisotropic quantum ground states 
in connection with generic phases of a doped Mott insula- 
tor have also been carried out£& It was suggested that as 
quantum fluctuation, induced by hole doping, increases, 
the Mott insulator turns into a smectic — stripe — phase, 
and a further increase of quantum fluctuation will turn a 
smectic to a nematic phase A£ The electronic smectic and 
nematic phases are characterized by broken translational 
(in one direction) and rotational symmetries, in analogy 
with classical liquid crystals^ 

The electronic nematic phase which breaks the point- 
group symmetry on a square lattice was studied in the 
extended Hubbard model. The instability of Fermi 
liquid towards possible ordered phases including the 
nematic — called the Pomeranchuk instabilitjM* — were in- 
vestigated, and it was shown that the nematic phase 
is the leading instability within some range of the pa- 
rameter spaceiii*i2ii^ Within the weak coupling theory, 
the effective Hamiltonian of the quadrupole-quadrupolc 
density interaction for the nematic phase was also con- 
structed, where the expectation value of the quadrupolc 
density is the order parameter of the nematic phased* 
The study of the nematic phase and possible supercon- 
ductivity in the continuum model within the weak cou- 
pling approach showed the non-Fermi liquid behavior in 
the single particle scattering ratei^ and an exotic super- 
conducting pairing symmetry^ via the coupling to the 
collective modes. Possible probes of the nematic phase 
were also discusscdiiSiii 

The nature of the quantum phase transition between 
the nematic phase and isotropic Fermi liquid as a function 
of chemical potential, for a particular set of parameters 
(interaction strength F2 and next nearest neighbor hop- 



ping t'), was recently studied numerically^ This compu- 
tation showed that the transition is generically first order. 
On the other hand, Ref . Il9l found anisotropic non-Fermi- 
liquid behavior at the isotropic-ncmatic quantum criti- 
cal point — second order transition — in the presence of a 
lattice. Therefore it is important to investigate whether 
the nature of isotropic-nematic transition depends on the 
type of electron dispersion, the interaction strength, or 
the temperature. 

In this paper, we investigate the behavior of the free 
energy and density analytically, as a function of chemical 
potential and interaction strength. We also extend the 
analysis to finite temperature and quasi-2D electron dis- 
persion. At zero temperature, we show that the nematic 
transition is first order as a function of the interaction 
strength and chemical potential for all typical 2D elec- 
tronic dispersions since they possess van Hove singular- 
ities. The transition changes to a continuous one at a 
finite temperature, but is not strongly affected by small 
dispersion in the third direction. 

We also find that the transition takes place for arbi- 
trarily small attractive interaction at the van Hove band 
filling. In the absence of interaction, at the van Hove 
filling, the Fermi surface changes topology from electron- 
to hole-like. This transition, first studied by Lifshita22*2I, 
causes singularities in thermodynamic quantities such as 
compressibility, due to van Hove singularities, but is not 
accompanied by any broken symmetry. This transition 
was recently revisited in Ref. |22j as an example of phase 
transition between different quantum orders which are 
not classified by broken symmetry. Our findings show 
that the Lifshitz transition docs not take place because 
the van Hove singularity is avoided due to a sudden 
change of the Fermi surface topology. This has implica- 
tions for earlier studies of broken symmetry states, such 
as density waves. These studie3ii*ii2i242Si2& were based 
on the existence of the van Hove singularity, which lead 
to a divergence of the relevant susceptibility indicating a 
transition to an ordered phase. Our results suggest that 
these ordered phases may be preempted by the first order 
transition into the nematic phase. 

The paper is organized as follows. We describe the 
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effective model Hamiltonian for the nematic in section 
im The mean field analysis at zero temperature and finite 
temperature is given in scction lTTTl A possible connection 
to other competing instabilities is discussed in section llVl 
We also provide the summary of our findings in the last 
section, IIVI 



II. MODEL 
A. Hamiltonian 

Our choice of model Hamiltonian is largely motivated 
by symmetry considerations and by the philosophy of 
the weak-coupling BCS theory. The interaction is cho- 
sen with foresight toward the mean-field analysis and 
with inspiration from classical liquid crystal theory. In 
a liquid crystal, each rod-like molecule defines a direc- 
tion in space, and the order parameter for the nematic 
phase is equally sensitive to their alignment as well as 
anti- alignment. It is a quadrupole (second order sym- 
metric traceless) tensor built from these directions. In 
two dimensions, it changes sign under a rotation by 90° 
and is invariant under a rotation by 180°. For an elec- 
tron gas, we can construct a similar order parameter 
from the momenta of electrons, the quadrupole density 
Qij = PiPj — ^P 2 ^ij- The interaction between quadrupole 
densities (suitably quantized and discretized) is made at- 
tractive to favor the alignment or anti-alignment of elec- 
tron momenta, i.e. formation of the nematic phase. 

The precise Hamiltonian under consideration is written 

as 



# = 5>k-M)4 c k + ^2(q)Tr QH<l)Q(<l) , (1) 



where £k is the single-particle dispersion, F 2 (q) is the 
inter-electron interaction strength, and Q(q) is the the 
lattice analog of the quadrupole density tensor con- 
structed in two dimensional square lattice from particle 
momenta given by 
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(2) 



The dispersion is that of a next-nearest-ncighbor tight 
binding model on a square lattice £k = — 2i(cos k x + 
cos k y ) — At' cos k x cos k y . While F 2 (q) is an arbitrary 
short ranged interaction, for example of the (two dimen- 
sional) Yukawa form, 



^(q) = - 



2 (1 + K 2 q 2 ) 3 / 2 ' 



(3) 



We shall only be interested in the strength of the inter- 
action F 2 , which is assumed to be positive. 




A/2t 



FIG. 1: Shapes of the Fermi surface (FS) for different values 
of /i and A. Lines (fJ,/2) 2 = A 2 mark the van Hove singu- 
larities in the DOS which correspond to the FS touching the 
boundary of the Brillouin zone. 



III. MEAN-FIELD ANALYSIS 

To decouple the quartic interaction in Eq. JQ) , we con- 
centrate on the q = interactions (which is equivalent 
to letting F 2 (q) — > — ^F 2 <5 q; o as n — ► 00 in Eq. Q) and 
define two order parameters, 

A = F 2 {Q xx (0)) and A' = F 2 (Q xy (0)), (4) 

where the expectation value (c^Ck) is replaced by the 
Fermi distribution ni?(ek — m); and £ k is the the renor- 
malized single-particle dispersion relation 

-2i(cos k x + cos k y ) — At' cos k x cos k y 

+ A(cos k x — cos k y ) — A' sin k x sin k y . (5) 

A finite t' does not qualitatively change the physics of 
the model. For simplicity, the value t' = was used in 
the T — calculations, but a nonzero t' was used at finite 
temperature. 

A non-zero expectation value of Qij indicates that the 
rotational (lattice point-group) symmetry has been bro- 
ken and a preferred direction for electron momenta has 
been selected. In particular, a non-zero expectation value 
of Q xy indicates that this direction is not parallel to ei- 
ther of the crystal axes. The profile of both the inter- 
action and the bare dispersion favor alignment along the 
crystal axes, as long as the coupling constants for diag- 
onal Q xx and off-diagonal Q xy elements are the same, 
F 2 X = F^ y . Hence we expect A' to vanish. This obser- 
vation has been confirmed by numerical calculations^ 

In general, unlike in the continuum modelia with full 



rotational symmetry, the coupling constants F^ x and 



F% v can be different. In such cases, a finite A' is possible. 
In fact broken symmetry states with oblique alignment, 
such as diagonal stripes, have been discussed in Ref. 

A. Free energy 

The resulting mean-field grand-canonical free energy 
density is given by 

F^,A) = y^+F ^,A), (6) 

where Fq is 

F (^A) =-~| deV(e)\n(l + e-^^), (7) 

with T>(e) being the density of states (DOS) of the single- 
particle dispersion (0. Its exact expression is (see Ap- 
pendix for details and definitions) 

(8) 

where A = l/(2tn 2 ). For |e| > 4t, V(e) vanishes. Note 
the logarithmic van Hove (vH) singularities in the DOS 
at (e/2) 2 = A 2 . These singularities occur when the con- 
stant energy contour (Fermi surface) meets the boundary 
of the first Brillouin zone as shown in Fig.^ Due to these 
remnants of the vH singularity of the non-interacting sys- 
tem, the free energy F(fi, A) will exhibit non- analytic be- 
havior at (/i/2) 2 = A 2 . This behavior is to be contrasted 
with the free energy in the presence of density wave order 
parameters, where all singularities are smoothed. Ulti- 
mately, persistence of non-analyticities in F(/x, A) leads 
to the first order isotropic-nematic quantum phase tran- 
sition. 



B. Zero temperature 
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fi/2t = n/2t = 0.089 n/2t = 0.109 




n/2t = 0.118 n/2t = 0.125 fj,/2t = 0.148 




FIG. 2: Plots of free energy as a function of A for different 
values of the chemical potential (F2N0 = 0.1) centered at 
A = 0; (b) /j, — (d) n — /i c (see text). 
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First we analyze the free energy density in the limit of 
zero temperature (/3 — > 00). In this limit 



F (^A) = (E-im), 



(9) 



where E and n are the energy and particle density per 
unit cell. 

The energy integral can be evaluated exactly, but the 
density integral cannot. The combined expression for the 
free energy (for A, # < 2t, and neglecting terms indepen- 
dent of A) is 
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FIG. 3: Phase diagram in the fj.-F2 parameter space. Solid 
curve, given by /j, c il l 2> . is the line of the first order transition. 



where for brevity all energy quantities are in units of 
2i. Any results extracted from this expression are valid 
up to quadratic order in fj, and A. The details of the 
calculations are outlined in the Appendix. 

The free energy for different chemical potentials is plot- 
ted as a function of A in Fig. [21 As the chemical potential 
decreases, the free energy develops local minima at finite 
A (the nematic phase), which then become the global 
minima for < /i c . It is clear that the transition be- 
tween isotropic and nematic phases is first order. 

From Eq. Q1U|). one can also show that A = is al- 
ways an extremum of the free energy. However, no local 
minimum of the free energy can be located in the region 
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FIG. 4: Equilibrium value of the order parameter A. The 
shaded region indicates negative curvature of the free energy. 
Lines A 2 = ) 2 are marked by dot-dashed lines. 



FIG. 5: The electron density per unit cell. The jump in the 
density is a signature of a first order phase transition. 



;f) 2 -A 2 | < (If) 2 , where 



H*/2t = 1.08e 4i? 2« 



fill 



this is the limit of metastability of the isotropic phase 
(A = 0) which becomes unstable for < /x*. Hence, A 
must have a finite equilibrium value in this region (ne- 
matic phase). The phase transition actually takes place 
slightly outside this region at = fi c , 



Hc/2t = 1.39e 4F 2 N o . 



(12) 



We stress again that \x c > yu» , which indicates that the nc- 
matic transition takes precedence over the Pomeranchuk 
instability (divergence of susceptibility). In Fig. we 
show the line of the first order phase transition in the 
parameter space of interaction strength, F%, and chemi- 
cal potential. 

The nontrivial local minima are located at 



±A/2t = 1.47e" 



- 1.36e« ( — 
\2t 



(13) 



Eq. ^|is valid for < /i c . The Fig. 0|shows the behav- 
ior of the order parameter. The order parameter jump 
and the width of the nematic window decrease exponen- 
tially as the coupling strength F 2 goes to 0. The unstable 
region |A 2 - (|) 2 | < (^) 2 is shaded in Fig. Q] A local 
minimum of the free energy can only be found outside 
this area, which explains the discontinuous character of 
the phase transition, since the equilibrium value of the 
order parameter must jump to avoid it. The dot-dashed 
lines, (^) 2 = A 2 , are where the change in topology of 
the Fermi surface takes place, the Lifshitz transition. As 
shown in Fig.0J these lines are embedded in the unstable 
region. Therefore, our results indicate that, in the pres- 
ence of interaction, the Lifshitz transition is not realized 
due to the formation of the nematic phase. 
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FIG. 6: The finite temperature phase diagram for the 
isotropic-nematic transition. The first order transition per- 
sists at low temperatures. However, the metastable region 
gets smaller and smaller with increasing temperature. Finally, 
around T/2t = 0.1, the phase transition becomes continuous. 



The density as a function of fj, is shown in Fig. No- 
tice that the dotted line is the density in the absence 
of nematic order. Its derivative at (i = is singular, 
which signals the Lifshitz transition. However, once the 
nematic order sets in, the density at half filling becomes 
smooth, but shows a discontinuity at \x c . 
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C. Finite temperature 

To investigate the robustness of the first order isotro- 
pic-nematic transition, we study the transition at finite 
temperature and with a finite inter-plane hopping term 
t z . These calculations are performed numerically using 
the same technique as in Ref.Hl The first order transi- 
tion is robust against a small t z = O.lt term, but changes 
to a second order one at a finite temperature. 

The phase diagram for finite temperatures is shown in 
Fig. El for a fixed F2N0 = O.fl and a finite next nearest 
neighbor hopping t 1 (t 1 = —OAt). The full width of the 
nematic window is from /i/2i = —0.94 to —0.57^ It is 
worth noting that the first order transition does not alter 
qualitatively with a finite t' . 

A negative (positive) t' shifts the window of the ne- 
matic phase to the hole- (electron-) doped side. At low 
temperatures the transition is still first order surrounded 
by a metastable region (which indicates the presence of 
unstable local minima). However, at about T/2t = O.f, 
the metastable region disappears and the transition be- 
comes continuous. 

The discontinuity in the isotropic-nematic transition 
at zero temperature can be traced to the presence of the 
lattice, which dictates the form of the dispersion relation 
and the presence of the van Hove singularity, and the 
sharpness of Fermi distribution. With increasing tem- 
perature, the thermal fluctuation will smear the sharp- 
ness of the Fermi distribution function, which results in 
a smaller jump in the order parameter and finally in a 
continuous transition. This expectation is confirmed by 
the numerics presented here. 



the understanding of the effective interaction for the ne- 
matic phase from the microscopic Hamiltonian is still 
missing, it suggests that the effective nematic interaction 
is hidden in the extended Hubbard model. 

Our discovery of a strong tendency toward the nematic 
phase near van Hove filling suggests that the two patch 
model should be revisited. Since the nematic order oc- 
curs for extremely small interaction near van Hove fill- 
ing, the existence of saddle points itself should be ad- 
dressed. Since the Fermi surface topology "suddenly" 
changes from closed to open as indicated in our result, 
we speculate that the instability towards competing or- 
ders such as antiferromagnetic and charge density wave 
(which are sensitive to the topology of the Fermi surface 
and enhanced by van Hove singularity) would be sup- 
pressed by the formation of the nematic phase. 

In summary, we have investigated a model Hamilto- 
nian exhibiting the nematic phase. At zero temperature 
the isotropic-nematic transition takes place for arbitrarily 
small coupling at van Hove band filling. Away from the 
van Hove filling, a finite minimum interaction is required 
to stabilize the nematic phase. The phase transition is 
first order as a function of interaction strength (chemi- 
cal potential) for a fixed chemical potential (interaction 
strength) as shown in Fig.|3J The strong tendency toward 
the nematic phase for an arbitrary small interaction at 
van Hove filling suggests that the Lifshitz transition is 
suppressed in the presence of interactions. At a finite 
temperature the transition becomes second order, while 
it remains first order for a quasi-2D dispersion. 
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IV. DISCUSSION AND SUMMARY 

There has been a number of one-loop renormalization 
group studies for the Hubbard model taking into account 
the Fermi surface only at the saddle points, (±7r, 0) and 
(0, ±7r), namely the two patch model at van Hove filling. 
These studies showed that there are antiferromagnetic 
and d-wave pairing instabilities in Hubbard modeliSi^SS 
Recently the two patch model was revisited, and trun- 
cation of the Fermi surface near the saddle points was 
suggested^ 

On the other hand, the instability of the Fermi li- 
quid towards the formation of the nematic phase — 
Pomeranchuk instability — with other competing orders 
were also recently investigated using different methods 
in the extended Hubbard modeliiiiSiiiiSii and the t-J 
models In Ref.Jn], the authors noticed that the nematic 
instability is driven by the attractive (repulsive) interac- 
tion between electrons in the same (different) patches via 
forward scatterings in the Hubbard model. This find- 
ing is consistent with our effective Hamiltonian, where 
— F 2 (cos k x — cos fe^) (cos k' x —cos k' y ) suggests an attractive 
interaction between electrons near (±7T, 0) and (±7r, 0), 
and a repulsive one between (±7r, 0) and (0, ±tt). While 
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APPENDIX A: DERIVATION OF FREE ENERGY 

The Density of states is defined and evaluated as 

c ( £ )4E*( £ - e *) 

k 

=N ^{v^w K { l -wi^w)}' (A1) 

where K(m) is the Complete Elliptic Integral of the First 
Kind^S and N = l/(27r 2 t). The function K{m) has a 
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logarithmic singularity at m = 1. At zero temperature, 
the free energy density is (cf. Eq. ©) 



1 A 2 

F = —— + (E - fxn). 



(A2) 



Here E is the energy density 

E= [ deeV(e) 
J— it 

= -(4<) 2 iVoRc| v /l-(ii) 2 i?(^l 



a 2 - m 



(A3) 



where E(m) is the Complete Elliptic Integral of the Sec- 
ond Kindi^S And n is the electron density 



-41 

(2tN ) 



(A-^)lnl 

V V 4 



A 



O -» -m) 



(A4) 



This expression is a leading order expansion in (A ± 
Here, for brevity, fj, and A are in units of It. 



Expanding the energy density to the same order, the 
combined free energy (neglecting terms independent of 
A) is given in Eq. (|10|l . Equating to zero the first deriva- 
tive of the free energy with respect to A gives an equa- 
tion for its local extrema. Nontrivial minima are easily 
found at (j, = 0. In Eq. (|13f) location of these minima is 
given to quadratic order in /i. Using this expression we 
can find the chemical potential \i c (Eq. (|12[) ~) at which 
they become global minima. Equating to zero the sec- 
ond derivative of the free energy with respect to A, we 
find the chemical potential ^* (Eq. Hill) ) at which the 
cxtremum at A = changes from a local minimum to a 
maximum. 
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